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S ■ Abstract 

The annihilation of electron-positron pairs around one second after the 
big bang distorts the Fermi-Dirac spectrum of neutrino energies. We 
determine the distortions assuming neutrino mixing with an inverted 
neutrino-mass hierarchy. Nonequilibrium thermodynamics, the Boltz- 
mann equation, and numerical integration are used to achieve the re- 
sults. The various types of neutrino behavior are established as a func- 
tion of masses and mixing angles. 
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Around a tenth of a second after the big bang, the early Universe was predomi- 
nantly a hot lepton-photon gas in thermal equilibrium Shortly thereafter, when 
the temperature of the Universe fell to about one MeV, the electrons and positrons an- 
nihilated into photons. The transfer of the lepton energy to the photon gas made the 
photon temperature slightly higher than if the lepton annihilation had not occurred. 
At about the same time, the neutrinos decoupled thermally because the frequency of 
hard-scattering processes decreased due to the lower neutrino energies. 

Since electron neutrinos couple more strongly to electrons and positrons than 
muon or tau neutrinos, the electron neutrinos received more energy from the electron- 
positron annihilations. In addition, higher-energy neutrinos decoupled later than 
lower-energy ones. As a result, distortions appeared in the neutrino Fermi-Dirac 
distributions. The timing of the neutrino decoupling is critical for the distortions to 
arise. If the neutrinos had decoupled much earlier or much later than the electron- 
positron annihilation, then little or no distortion of thermodynamic distributions 
would have arisen. One consequence of the timing coincidence is that an excess of v e 
over and v T formed around this time. The distortions have been analyzed in the 
case of no neutrino mixing in refs. [fj. Q . 

The issue of how the standard model of the early Universe is affected by the 
presence of various types of neutrino and of neutrino masses and mixings has been 



investigated by several authors during the past few years P-[|16fl. In the present 
work, we assume that only electron and muon neutrinos have significant mixing. 
Then, on average no neutrino oscillations occur in the early Universe before the v e 
excess develops, because until then for every v e that oscillates into a there is a 
that oscillates into a v e . However, neutrino mixing can affect the generation of the 



spectral distortions [fLq| . Also, after the neutrinos decouple and the distortions are 
established, energies and densities continue to fall as the Universe expands and so the 



relevance of neutrino mixing increases JTT|, [14| . 

The neutrino mixing can be characterized by the vacuum neutrino-mixing angle 9 
and the difference A = m\ — m\ between the squared masses of the mass eigenstates. 
The parameter region A < corresponds to an inverted neutrino-mass hierarchy, in 
which the heavier mass eigenstate has a larger component in the electron neutrino 
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than in the muon neutrino. For A < 0, the flavor evolution of neutrinos in the early 
Universe exhibits a rich spectrum of features [ I2[ . Several recent papers have been 
based on the possibility of an inverted neutrino-mass hierarchy |T7|, IB|, [L9| . 

For |A| < 1CT 6 eV 2 , neutrino mixing effects are too small to affect significantly 
the production of the v e excess, so neutrino oscillations are relevant only after produc- 
tion. It is therefore a good approximation to assume the production phase, t < 0.33 
seconds, generates the known zero-mixing distortions and to adopt the results as 
input for the oscillation phase [JTT], 0, [l4| . This procedure avoids the difficulties of 
simultaneously treating the nonequilibrium neutrino-excess production and the neu- 
trino oscillations. 

In the region 10~ 6 eV 2 < |A| < 10~ 5 eV 2 , oscillations play some role in the dis- 
tortion of neutrino distributions, while for |A| > 1CU 5 eV 2 oscillations definitely must 
be incorporated in the excess production process. For A > 1CU 6 eV 2 , another ap- 
proximation is available to bypass the complexities of the full production-oscillation 
process: the appropriate initial conditions for the pure-oscillation phase can be es- 
timated by projecting the zero-mixing neutrino excess onto approximate nonlinear 
mass-eigenstate (ANME) configurations |T2|, [TJJ]. It has recently been demonstrated 
that ANME configurations are indeed the correct initial conditions for the oscillation 
phase for the case of positive A > 10~ 6 eV 2 |TB . 



To date, no analysis of neutrino oscillations in the early Universe has been per- 
formed for negative A < — 1CT 6 eV 2 . The basic reason is that early- Universe neutrino 
oscillations are unstable when A is negative, and so at some point in time the av- 
erage neutrino flavor undergoes a dramatic change in behavior. This time is called 
the instability time in what follows. When — 1CT 6 eV 2 < A < 0, the instability time 
occurs after the production phase, and hence does not affect the neutrino excess. The 
method of ANME configurations can therefore be used [Q. However, for A < — 10~ 6 
eV 2 , the instability time occurs during the production phase, which makes it difficult 
to surmise appropriate initial conditions for the subsequent neutrino oscillations. 

The primary purpose of the present work is to analyze this remaining unexplored 
parameter region. Using the full apparatus for the analysis of combined oscillations 
and production, we have obtained results for the neutrino spectral distortions and 
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for neutrino oscillations for A < —10 6 eV 2 . A secondary purpose is to verify the 



validity of the ANME-based initial conditions assumed in ref. [12[] for the oscillation 
phase in the region — 10~ 6 eV 2 < A < 0. 

The lengthy equations governing the excess-^ e production in the presence of os- 
cillations may be found in section V of ref. [ IE| . They are based on nonequilibrium 
thermodynamics using the Boltzmann equation, and they incorporate hard-scattering 
production effects from interactions among photons, electrons, positrons, and neutri- 
nos, as well as neutrino forward scattering off neutrinos, electrons, and positrons. 
Neutrino-neutrino forward scattering is often the dominant neutrino-oscillation effect 
by many orders of magnitude. It renders the system of equations nonlinear, so nu- 
merical simulations are necessary. In the present work, we employ the notation and 
methods of ref. fL6| . In particular, v (E, t) = [y\v e + v^v^) is the sum of the excess of 
electron and muon neutrinos at energy E per canonical comoving volume per E/T, 
Vi(E,t) = [y\v e — u^i/p) is the difference, v^E^t) = v\v T is the excess of tau neutri- 
nos, while V2(E,t) = 2Re{y\v p ) and vz(E,t) = 21m{y\v p ) are the real and imaginary 
parts of the off-diagonal component of the density matrix for neutrinos with energy 
E. In units with h — c — 1, the canonical comoving volume is taken to be 1.0 MeV -3 
when the temperature T of the Universe is 1.5 MeV. 

For positive A, a single behavior characterizes neutrino flavor evolution [X?|: al- 
though the flavor may change over time, any evolution occurs smoothly. In contrast, 
when — 10 -6 eV 2 < A < three types of behavior appear [T2j depending on the 
values of 9 and A. Near sin 2 29 = 1, the evolution is smooth. For sin 2 29 small, self- 
maintained coherence occurs. This is a mode exhibiting some features reminiscent 
of solitons, in which a fraction of the neutrinos oscillate essentially in phase. It was 
originally observed in test simulations [[HJ and is now understood mathematically 
fHfl . The third behavior, called irregular in ref. [12], occurs for intermediate sin 2 29 
and relatively large |A|. 

Mathematical and physical insight about irregular behavior has recently been 



obtained [22]. Ignoring the expansion of the Universe but incorporating the alignment 
property of the neutrino vectors [[12], H| , analytical solutions for the behavior of the 
average neutrino flavor were found. These solutions are characterized by two modes. 
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When both have very small amplitude, smooth behavior occurs. When only one mode 
has sizeable amplitude, self-maintained coherence arises. Finally, when both modes 
are significantly excited, bimodal self-maintained coherence appears. This can mimic 
the irregular behavior seen in ref. [12]. It has been conjectured that, when the periods 
of the two modes are incomensurate, bimodal self-maintained coherence is irregular 
behavior [^2] . 

To verify that incorporating hard-scattering production effects leaves unaffected 
the results of ref. |12|], we performed full production-oscillation simulations for se- 
lected values of sin 2 29 and A in the region — 1CT 6 eV 2 < A < 0. These simulations 



qualitatively confirm the results of ref. [12] throughout the entire parameter region. 
They also agree quantitatively over all regions displaying self-maintained coherence 
or smooth evolution. For example, Figures 2 and 4 of ref. |12| are reproduced. In the 
irregular region, the numerical results are qualitatively the same, and the quantita- 
tive discrepancy can be understood in terms of bimodal coherence. It turns out that 
the effect of incorporating the true neutrino-excess production can excite the ampli- 
tudes of the two modes by slightly different amounts. Since the periods are then 
incommensurate, small differences can lead to sizeable numerical changes. Bimodal 
coherence exhibits features similar to chaotic behavior, in that it is very sensitive to 
small changes in parameters and initial conditions. 

To determine the behavior of a particular simulation, one can examine the time 
evolution of the components of average neutrino vector (v(t)) = ((t>i(t), t> 2 (£), t>3 (£))). 
However, the single or bimodal nature of coherent oscillations is often better displayed 
as the three-dimensional orbit of (v(t)). In the current work, bimodal oscillations 
were evident in almost all the simulations in the irregular region. Figure 1 shows an 
example. In a few of the simulations, the orbits appear irregular: it is impossible to 
decide if only two modes occur. It may be that some irregular behavior is not strictly 
bimodal coherence due to the expansion of the Universe or other effects not included 



in the analysis of ref. [22 . 

Figure 2 shows our results in the sin 2 29- A plane. The region labelled 'weak 



bimodal,' classified as 'residual self-maintained coherence' in ref. [[12[], is bimodal 
coherence with one mode having a much smaller amplitude than the other. It is almost 
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self-maintained coherence but exhibits irregular-type behavior over sufficiently long 
time scales. Simulation for A < — eV 2 is prohibitive in computer time, but we 
can determine the behavior by other arguments. The essential physics is unchanged 
from A = — 1CT 4 eV 2 to A = — 1 eV 2 , so irregular behavior should appear in this 
region. For A < — 1 eV 2 , vacuum oscillation effects dominate over the nonlinear 
forward scattering of neutrinos, and so vacuum behavior is expected. This region 
is labelled 'vacuum decoherence' in Figure 2. The production of excess neutrinos 
quickly decoheres and neutrinos become vacuum mass eigenstates. 

For positive A, the sum of the excess of electron and muon neutrinos is insensitive 



to the effects of neutrino oscillations [16) . The same is true of the tau-neutrino excess. 
For negative A, we find similar results. Hence, our new curves for (vo(t)) and (v^(t)) 
exhibit little dependence on A and sin 2 29. They are essentially identical to the 
dashed-dotted and dashed curves in Figure 3 of ref. [|16[]. Likewise, the final energy 
distributions for the sum of the excess of electron and muon neutrinos and for the tau- 
neutrino excess are insensitive to A and sin 2 29 and are given by the dashed-dotted 



and dashed curves in Figure 1 of ref. |L6 



In contrast, the remaining three variables (vi(t)), (t> 2 (t)), and (v 3 (t)) depend on 
A and sin 2 29. The difference (vi(t)) between the excesses of electron and muon 
neutrinos starts out at zero at t — 0.08 seconds when T = 3.0 MeV, rises as excess 
electron neutrinos are produced, and typically declines after the instability time as 
many electron neutrinos oscillate into muon neutrinos. At the instability time and 
thereafter, (v2{t)) and sometimes (vs(t)) become significantly different from zero. 
The sizeable conversion of electron neutrinos into muon neutrinos throughout the 
A < region is the result of a new resonance conversion mechanism |L2]| driven by 
the imaginary part of the off-diagonal component of the neutrino density matrix, 
which represents the difference between the third components of the average vectors 
for neutrinos and antineutrinos. 

Figure 3 displays (vi(t)), ^(t)) and (i>3(t)) for the case A = — 10~ 5 eV 2 and 
sin 2 29 = 10~ 3 , a point lying in the irregular region of Figure 2. Figure 4 displays the 
same three quantities for the case A = — 10~ 6 eV 2 and sin 2 29 = 10 -8 , which lies in 
the weak-bimodal region. Rapid fluctuations occur that are not visible on the scale of 
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Figure 4. In Figure 5, we plot the components over a short time interval at t ~ 0.94 
seconds. The corresponding orbit is shown in Figure 6. Since irregular behavior 
occurs throughout most of the region with A < — 10 -6 eV 2 and vectors undergo 
sizeable unpredictable fluctuations, it makes little sense to plot the final production 
profiles for v as a function of energy. 

For — 10~ 4 eV 2 < A, ANME configurations at early times are flavor eigenstates 
with an excess of electron neutrinos, so the first component of v(E,t) is the largest. 
At later times, ANME configurations undergo a 180° rotation in three-space. If the 
nonlinear neutrino- neutrino interactions were absent, then little conversion of electron 
neutrinos into muon neutrinos would occur when the rotation is non-adiabatic. How- 
ever, after the rotation occurs, the neutrino-neutrino interactions render the system 
unstable and a large conversion of electron neutrinos into muon neutrinos takes place, 



regardless of the rotation rate |12|, [T4]|. The instability time is expected to occur when 
the vacuum and CP-conserving electron-neutrino terms become equal [JT2] , |I~4}1 , i.e., 
when 



| A| cos(20)/(2£) « \V CP+ \ = 2V2G F E (p e - + p e - + p e+ + p e+ ) /M\ 



w 



where E is the average neutrino energy, Gp is Fermi coupling constant, M\y is the 
mass of the W boson, and p e ± and p e ± are, respectively, the energy densities and 
pressures of the background positrons and electrons. 

When 9 is small so that cos 29 « 1, we find that Eq. (p]) predicts instability 
times of 1.1, 0.5, and 0.23 seconds for A = -10~ 7 eV 2 , -10" 6 eV 2 and -10~ 5 eV 2 , 
respectively. In our numerical simulations, the corresponding instability times are 
found to be ~ 1.0, ~ 0.45 and ~ 0.2 seconds, in good agreement with the theoretical 
expectations. 

Similarly, when sin 2 29 is near one we again find good agreement. For example, 
when sin 2 29 = 0.49, the instability times occur in simulations at ~ 0.55, ~ 0.25 and 
~ 0.12 seconds for A = -10~ 6 eV 2 , -10" 5 eV 2 and -10~ 4 eV 2 , respectively. The 
corresponding values from Eq. (|l|) are 0.52, 0.24 and 0.11 seconds. 

For A < — 10 -4 eV 2 , the instability time as determined by Eq. ([I]) occurs before 
production begins. Hence, the instability appears as soon as the electron-neutrino 
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excess becomes sufficiently large to render the neutrino-neutrino interactions impor- 
tant. This happens soon after production begins. The instability disrupts the system 
and leads to the irregular behavior seen in this region in Figure 2. 

Other features of the neutrino oscillations are similar to the positive- A case [14]. In 



particular, the CP-suppression mechanism is still operative. In the regions of coherent 
oscillations, the CP asymmetry of the excess can sometimes reach the 1% level, but 
this is insufficient to affect big-bang nucleosynthesis at a level observable with current 
technology. The evolutions of the antineutrino excesses are almost identical to those 
of the corresponding neutrinos, so we have not displayed them. 

To summarize, we have obtained the complete phase diagram for neutrino oscil- 
lations for the case of an inverted neutrino-mass hierarchy. Five different behaviors 
are observed: smooth evolution, self-maintained coherence, bimodal coherence, irreg- 
ular behavior, and vacuum decoherence. Throughout most of the large-|A| region, 
irregular coherent oscillations are seen. 
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Figure Captions 



Figure 1: An orbit in the irregular region for sin 2 29 = 1CT 7 , A = — 10 _5 eV 2 near 
t = 0.42 seconds. The first mode involves motion around a circle. The second mode 
rotates the circle. Only a few periods of the first mode are included. When the orbit 
is displayed over longer time intervals, the bimodal coherence produces a figure which 
looks like a ball of wool. 

Figure 2: The phase diagram for neutrino flavor oscillations for an inverted neutrino- 
mass hierarchy. 

Figure 3: Time development of the nonequilibrium distortions for the case A = — 10~ 5 
eV 2 and sin 2 29 = 10~ 3 . (a) The quantity (vi(t)) is shown as a function of time, (b) 
The quantity (^(i)) is shown as a function of time. Note the change of vertical scale 
relative to (a), (c) The quantity (v 3 (t)) is shown as a function of time. Note the 
change of vertical scale relative to (a). 

Figure 4: Time development of the nonequilibrium distortion for the case A = — 10~ 6 
eV 2 and sin 2 29 = 10~ 8 . (a) The quantity (vi(t)) is shown as a function of time, (b) 
The quantity (^(i)) is shown as a function of time, (c) The quantity (v 3 (t)) is shown 
as a function of time. 

Figure 5: Time development of the nonequilibrium distribution for the case in Figure 
4, with A = — 1CT 6 eV 2 and sin 2 26 l = 1CT 8 . The components have been scaled by a 
factor of about 250 for clarity. 

Figure 6: Orbit for the case A = — 10~ 6 eV 2 and sin 2 29 = 10 -8 . The three quantities 
displayed in Figure 5 are combined into a vector and its trajectory is plotted in 
three-dimensional space. 
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Figure 4a: first-component distortion 
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This figure "figure2.gif" is available in "gif" format from: 
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Figure 3b: second-component distortion 
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Figure 4b: second-component distortion 
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Figure 3c: third-component distortion 
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Figure 4c: third-component distortion 
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Figure 3a: first-component distortion 
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